hmax=9;
h=-hmax:hmax;
m=-hmax:hmax;
for i=1:(2*hmax+1)
    m(i)=mut_ind(h(i));
end
plot(h,m)